ON ABSORBING BOUNDARY CONDITIONS FOR LINEARIZED EULER EQUATIONS 

BY A PERFECTLY MATCHED LAYER 


Fang Q. Hu 

Department of Mathematics and Statistics, Old Dominion University 

Norfolk, VA 23529 


ABSTRACT 

Recently, Berenger introduced a Perfectly Matched Layer (PML) technique for absorbing elec- 
tromagnetic waves. In the present paper, a perfectly matched layer is proposed for absorbing 
out-going two-dimensional waves in a uniform mean flow, governed by linearized Euler equations. 
It is well known that the linearized Euler equations support acoustic waves, which travel with the 
speed of sound relative to the mean flow, and vorticity and entropy waves, which travel with the 
mean flow. The PML equations to be used at a region adjacent to the artificial boundary for 
absorbing these linear waves are defined. Plane wave solutions to the PML equations are devel- 
oped and wave propagation and absorption properties are given. It is shown that the theoretical 
reflection coefficients at an interface between the Euler and PML domains are zero, independent 
of the angle of incidence and frequency of the waves. As such, the present study points out a 
possible alternative approach for absorbing out-going waves of the Euler equations with little or 
no reflection in computation. Numerical examples that demonstrate the validity of the proposed 
PML equations are also presented. 
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1. INTRODUCTION 


Numerical boundary conditions have always been an important issue in computational fluid dy- 
namics. In problems involving wave radiation and convection, such as in computational acoustics, 
the physical domains are often necessarily truncated due to the limitation of a finite computational 
domain. Thus, at these artificial boundaries, numerical non-reflecting or absorbing boundary con- 
ditions are needed so that the out-going waves are not reflected. Various computational techniques 
have been developed in the past to minimize the reflection of out-going waves. They include the 
out-flow boundary conditions based on the characteristics of the Euler equations [1, 2] and radiation 
boundary conditions based on the far field asymptotic solutions [3, 4, 5, 6]. A recent review can be 
found in reference [7]. In addition, a buffer zone technique has been developed in which the mean 
flow is altered gradually to be supersonic in a buffer region adjacent to the artificial boundary [8, 
9], 

Recently, Berenger introduced a Perfectly Matched Layer (PML) technique for absorbing elec- 
tromagnetic waves in a finite difference method of solving the Maxwell equations [10]. In this 
approach, a PML medium of certain depth is introduced in a region adjacent to the artificial 
boundary of a computational domain. The equations for the PML medium are designed such that 
the out-going electromagnetic waves are absorbed by the layer with no reflection (theoretically). 
The novelty of the technique lies in the way the PML equations are constructed. Numerical studies 
have shown that the PML technique dramatically reduces the reflection of out-going waves when 
compared to several previous methods in electromagnetics [10, 11]. 

In this paper, a perfectly matched layer is proposed for absorbing out-going waves of the 
linearized Euler equations for two-dimensional problems. It is well known that the Euler equations 
linearized around a uniform mean flow support acoustic waves as well as vorticity and entropy waves. 
The acoustic waves travel with the speed of sound relative to the mean flow, while the vorticity 
and entropy waves convect downstream with the mean flow. To apply the PML technique, the 
computational domain is divided into the interior domain, where the Euler equations are used, 
and PML domains adjacent to artificial boundaries where the proposed PML equations are to be 
used (Figure 1). The purpose of the PML domain is to absorb the out-going waves. It will be 
shown that the theoretical reflection coefficients for incident linear waves at an interface between 
the interior domain and a PML domain are zero and the amplitudes of waves that enter the PML 
domain decrease exponentially. It will also be shown that it is true independent of the angle of 
incidence and frequency of the waves. 

In the next section, equations for the PML domains are defined. In section 3, plane wave 
solutions of the PML equations are developed and wave propagation and absorption properties are 
given. Then, the theoretical reflection and transmission coefficients at interfaces are calculated in 
section 4. Numerical examples that demonstrate the validity of the proposed PML equations are 
provided in section 5. Section 6 contains the conclusions. 
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2. PERFECTLY MATCHED LAYER EQUATIONS 

We consider the two-dimensional linearized Euler equations with a uniform mean flow : 
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in which u and v are the velocity components in the x and y directions respectively, p is the pressure, 
and p is the density. The velocities have been non-dimensionalized by the speed of sound a, the 
density by p and the pressure by pa 2 , where p is the mean density. For simplicity, a mean flow of 
Mach number M in the direction of the avaxis has been assumed. Situations of a mean flow not 


aligned with the avaxis are discussed in the Appendix. It is also assumed that the mean flow is 
subsonic, i.e. , M < 1. 

Equations (1.1)- (1.4) support acoustic waves, which travel with the speed of sound relative to 
the mean flow, and vorticity and entropy waves, which travel with the mean flow. Our aim is to 
define a perfectly matched layer to be used at a region adjacent to the artificial boundary which 
absorbs the out-going acoustic, vorticity and entropy waves with little or no reflection. 

Different from [10], here we split u, v, p and p in equations (1.1)-(1.4) into sub-components 
Ui, u 2 , Vi, v 2 , pi, p 2 and pi, p 2 . We define the following equations for the perfectly matched layer 
(PML) : 
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In the above, a x and a y have been introduced for the absorption of waves in the layer. They 
will be called absorption coefficients in this paper and are assumed to be greater than or equal 
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to zero. It is to be noted that, when a x = a y = 0, equations (2.1)-(2.8) can be reduced to the 
Euler equations (1.1)-(1.4) with u = Ui + u 2 , v = tq + v 2 , p = Pi + p 2 and P = Pi + p 2 . Thus the 
Euler equations are a special case of the PML equations. Moreover, the spatial derivatives involve 
only the total u, v, p and p which are assumed to be continuous at interfaces. We also note that 
two kinds of interfaces are created, namely, the interfaces between the interior domain and a PML 
domain and those between two PML domains, as shown in Figure 2. The former, of course, can be 
regarded as a special case of the later. In the next two sections, we show the wave propagation and 
absorption properties within a perfectly matched layer defined above and calculate the reflection 
and transmission coefficients at an interface between two PML domains. 


3. PLANE WAVES IN A PERFECTLY MATCHED LAYER 

Let a plane wave in the PML domain be expressed as 

[u u u 2 ,v u v 2 ,p u p 2 , p u p 2 \ = [n 10 , u 20 , v 10 , v 20 ,p l0 ,p 20 , p l0 , (3) 


in which a subscript 0 has been used to denote the amplitudes of the components. By substituting 
(3) into (2.1)-(2.8), we get 


(lo + ia x )u 10 = k x (p lQ + p 2o ) (4.1) 

(u + ia x )u 20 = k x M(u 10 + u 20 ) (4.2) 

(to + i(T y )v 10 = k y (p l0 + p 2o ) (4.3) 

(u + ia x )v 20 = k x M(v 10 + v 20 ) (4.4) 

(lo + icr x )pi 0 = k x (u 10 + u 20 ) + k x M(pi 0 + p 2o ) (4.5) 

(lo + ia y )p 2o = k y (v 10 + v 20 ) (4.6) 

(ll> + io'x)Pio = k x (ui 0 + u 20 ) + k x M(pi 0 + p 2o ) (4-7) 

(lo + ia y )p 2o = k y (v 10 + v 20 ) (4.8) 


It will be shown below that equations (4.1)-(4.8) support acoustic waves, when uj — k x M + ia x 7 ^ 
0, and vorticity and entropy waves, when u — k x M + ia x = 0. 


3.1. Acoustic Waves 

When lo — k x M + ia x 7 ^ 0, it is easy to find that the amplitudes of the components in (3) can 
be expressed in terms of n 10 and n 10 as follows : 
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In addition, from (4.1) and (4.3), we have the following relation : 


uj + ia x u 10 _ k x 
uj + io y v 10 k y ’ 


Now, by substituting (5.3)-(5.4) into (4.1) and (4.3), we get 


(uj + ia x )u 10 
(uj + ia y )v 1 o 


k 2 x (uj + icr x )(uj + i(T y )u 10 + k x k y (uj + ia x fvi 0 
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(7.1) 
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respectively. For (7.1) and (7.2) to have non-trivial solutions for n 10 and n 10 , it is found that the 
following dispersion relation for k x , k y and uj has to hold: 


(lj - k x M + ia x f(uj + ia y f - k 2 x (uj + ia y f - k 2 y (uj + ia x f = 0. 


(8) 


However, it has been found more convenient to express k x and k y in terms of n 10 and n 10 . On 
eliminating k y and k x in the numerators of (7.1) and (7.2), respectively, using equation (6), we 
obtain 
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(9.1) 

(9.2) 


The positive and negative signs indicate the direction of wave propagation. The positive sign will 
be taken in the discussions followed. For convenience, we express n 10 and n 10 as 


Ui o = A cos <f>, (10-1) 

v 10 = A sin <j>, (10.2) 

where A is a complex number and cj) is a real number. This is possible when n 10 and n 10 have the 
same phase, i.e. , when ii 10 /ti 10 is real. It will be shown in Appendix A that iiy/bo i s indeed real 
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for solutions of the Euler equations. Thus, it is sufficient to consider only solutions where w 10 /u 10 
is real. 

Substituting (10.1)-(10.2) into (9.1)-(9.2) and solving for k x and k y , we get 
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cos <f>, 


( 11 . 1 ) 


U + IGy 
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( 11 . 2 ) 


With these expressions, further simplifications can be made for the amplitudes in (5.1)-(5.6). As a 
result, the plane wave solution to (4.1)-(4.8) can be expressed as 
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It is easy to see that the above expression represents a wave propagating with the speed of 
sound (which is unity in the non-dimensional variables) relative to the mean flow in the direction 
making an angle cj) with respect to the avaxis [12]. This solution, thus, represents the acoustic wave 
in the PML domain. Furthermore, when a x or a y is not zero, the magnitude of the wave decreases 
exponentially as it propagates in the x or y direction respectively. 


3.2. Vorticity and Entropy Waves 

When lo — k x M + ia x = 0, it follows immediately that 


k x = 


v + io x 
M 


(13) 


Under this condition, it also follows that n 10 = 0, n 10 = 0 and p l0 +P 20 = 0 by (4.2), (4.4) and (4.1) 
or (4.3), respectively. Equations (4.1)-(4.8) now admit vorticity-like and entropy-like solutions. For 
convenience, we express n 20 , ^20 an( i Pi 0 as 


“20 = ~B sin -0, (14.1) 

v 20 = B cos -0, (14-2) 


Pi 0 — C. 


(14.3) 
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Then, k y for the vorticity wave is found to be 
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and the plane wave solution to (4.1)-(4.8) can be expressed as 
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We remark that the constants B and C independently represent the vorticity and the entropy 
waves. Since each component is of the form f( jp — /, y)e~^ x , the waves are convected with the 
mean flow at a speed that is equal to M . Furthermore, when a x is greater than zero, the magnitude 
is decreased exponentially in a;. It is worth pointing out that the vorticity waves do not contribute 
to the Ui an( i v i components in the PML equations. 


4. REFLECTION AND TRANSMISSION AT AN INTERFACE 
BETWEEN TWO PERFECTLY MATCHED LAYERS 

We now consider the wave reflection and transmission at an interface between two PML do- 
mains. This, of course, includes the interface between the interior domain and a PML domain. As 
in the PML technique for electromagnetic waves [10], the absorbing coefficients, a x and a y , will be 
chosen such that a y is the same across an interface normal to x and a x is the same normal to y. 
Since the Euler equations for the interior domain can be considered as PML equations with both 
absorption coefficients being zero, a y or a x will be consequently zero across an interface normal to 
x or y between an interior domain and a PML domain. This is as shown in Figure 2. In what 
follows we show that the reflection coefficient at an interface downstream normal to x is zero for 
incident acoustic, vorticity and entropy waves. Similar results can be established analogously for 
other interfaces. 

Let the interface be located downstream at x = 0 and the absorption coefficients be a xX and 
o y on one side and a x2 and a y on the other (figure 3). For a subsonic mean ffow, the possible 
reflection is an acoustic wave. Then, by the results of previous section, the incident, reflected and 
transmitted waves can be expressed as follows : 

(a) incident wave : 
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(b) reflected wave : 
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c) transmitted wave : 
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The angles of the acoustic waves are as indicated in Figure 3. 

At the interface, we impose the condition that ui+u 2 , v 1 + v 2 , P 1 +P 2 and P 1 AP 2 be continuous. 
Since this continuity is true for all values of y along the interface, it follows that the coefficients of 
y in the exponents of (16.1)-(16.3) must be the same for all the incident, reflected and transmitted 
waves. This yields 

sin fa sin fa 

1 — M cos fa 1 + iff cos fa ’ 
sin fa sin fa 

1 + iff cos fa 1 + iff cos fa ’ 
tan fa = tan fa. 

(Here, tan fa = sin<^/(l + M cos fa) has been assumed for the incident waves.) From above, it is 
found that 


cf> r 

4>t 

i’t 


2 tan 1 

4*i 1 

fa- 


1 - M 
1 + M 
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Furthermore, by the continuity of Ui + u 2 , rq + v 2 , Pi + p 2 and Pi + /° 2 5 we have 


(17.1) 

(17.2) 

(17.3) 


Ai(l + M cos fa) cos fa — A r ( 1 — M cos fa) cos fa — sin fa = A t ( 1 + M cos fa) cos fa — B t sin fa, 


Ai( 1 + M cos fa ) sin fa + A r ( 1 — M cos fa ) sin fa + cos fa = A t ( 1 + M cos fa ) sin fa + B t cos fa , 
Ai ( 1 + M cos fa ) + A r ( 1 — M cos fa ) = A t ( 1 + M cos fa ) , 

and 

Ai ( 1 + M cos fa ) + A r ( 1 — M cos <(y ) + Bi ' ^ 1 ~ Ci = A t ( 1 + M cos fa ) + B t ' ^ \~ Ct ■ 

Applying (17.1) and (17.3), the above equations can be rewritten as a linear system of four homo- 
geneous equations for A t — A t , A r , B t — B t and Ci — C t and the coefficient determinant is found to 
be 

2(1 + M cos fa)(l — M cos fa) cos ( ^ ^ < ^ r ^ cos (fa — ^ 1 < ^ r 


. 2 J 2 

It is straightforward to verify that the above is not zero for any angle of incidence and, thus, the 
only solution is 


A r = 0, 

(18.1) 

At = Ai, 

(18.2) 

Bt = Bi , 

(18.3) 

C t = Ci. 

(18.4) 
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Therefore, equations (17.1)-(17.3) and (18.1)-(18.4) demonstrate that at an interface between 
two PML domains downstream normal to the avaxis with absorption coefficients (a x i,a y ) and 
(cr X 2 i (J y) respectively, the reflection is null and the transmitted waves maintain the same direction 
and amplitude as the incident waves at the interface. This has been shown to be independent of 
the angle of incidence and frequency of the waves. 


5. NUMERICAL EXAMPLES 

Three numerical examples are presented in which the mean flow is in the direction of the x- 
axis or at a 45° angle with the avaxis and Gaussian shaped acoustic, vorticity and entropy pulses 
are initiated in the interior domain at / = 0. The first two examples have four open boundaries, 
while the third example has a solid wall at one boundary. The purpose of these examples is to 
demonstrate the validity of the proposed PML equations. 

The implementation of the PML equations in a finite difference scheme can be straightforward, 
since the spatial derivatives involve only the total u, v, p and p. A central difference scheme has 
been used for spatial discretization in the examples. In particular, a 4th-order 7-point explicit 
central difference scheme is adopted [6]. Time integration is carried out by a 4th-order Runge- 
Kutta scheme optimized for low-dissipation and low-dispersion errors [13]. Since a wide stencil has 
been used here, the absorption coefficients are varied gradually in the PML domain. Specifically, 
the following form has been used for a x and a y in the calculations presented : 


® 


d 

D 


P 


•> 


where D is the thickness of the PML domain and d is the distance from the interface with the 
interior domain. 

At the end of a PML domain, certain boundary conditions, such as a solid wall condition or 
other radiation boundary condition, can be applied. If this is the case, the wave is reflected and, 
when it re-enters the interior domain, the total absorption factor can be estimated as 


2 & rn -D cos <p 

/3 + 1 l_M 2 cos 2 4> 


for the acoustic wave in a PML domain normal to x [10]. The absorption factor for the vorticity 
and entropy waves is estimated as 

cr m D 

e ~ M(fi + 1 ) _ 

The parameters of the layer can thus be adjusted for desired absorption. 

In the first two examples, however, since the numerical solutions are decaying exponentially 
toward the edges of the computational domain, due to the presence of four PML domains, a periodic 
boundary condition is used at the end of the PML domains for computing the spatial derivatives. 
This eliminates the need of backward differences. Consequently, this will result in the out-going 
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wave traveling in the two opposite PML domains and the absorption factor for the acoustic wave 


is now 

2cr m D cos 4> 

g /3 + 1 1 ±M cos 4> 


where the + and — signs are for the right and left traveling waves respectively. We note that the 
absorption factor for acoustic waves varies with the wave angle. There will be, for instance, little 
absorption in a layer normal to x for waves with cj) close to ±7t/2. These, however, are absorbed 
by the layers normal to the y - axis [10]. Moreover, since the solution is treated as periodic in x and 
y , spectral filtering with FFT is employed. 

The initial conditions are adopted from the Benchmark Problems of Computational Aeroa- 
coustics [14] : 


p = e - ( ln2 ) — 


i ) J + (y-ya r 


O.le 


— 0 n 2) 


(x-x b ) J +(y-y b ) 


p = e -(ln2) <— +<«-«»> 


and 


u = 0.05(2/ - Vb)e (ln2) 


(x-x b ) 2 +(y-y b ) 2 


v = — 0.05(A — x b )e ( - ln2 - ) 


(x-x b ) 2 + (y-y b ) 2 


The above equations include an acoustic pulse centered at ( x a ,y a ) and a vorticity and an entropy 
pulse both centered at (x b ,y b ). 

In the first example, the mean flow is in the direction of the avaxis with M = 0.5. The 
computational domain is [—60,60] X [—60,60] where Ax = Ay = 1 and a PML domain of 10 grid 
points has been used around the four sides. Thus, the interior domain where the Euler equations 
are applied is [—50,50] X [—50,50]. The initial center of the acoustic pulse is ( x a ,y a ) = ( — 25,0) 
and that of the vorticity and entropy waves is (x b ,y b ) = (25,0). For the present calculations, 
a m = (3 = 2 has been taken. The contours of computed p and u are shown in Figure 4 and 5 for 
time t = 10, 40, 50 and 70. In Figure 6, contour plots of p within the PML domain are shown 
for the right boundary. Clearly, no visible reflection is seen in the contour plots as the pulses exit 
the boundaries. Comparisons with the exact solution are shown in Figure 7 for density along the 
line y = 0. Good agreement is observed. In the second example, the mean flow is in the direction 
making a 45° angle with the avaxis and M x = M y = 0.5 sin(7r/4). The initial locations of the pulses 
are ( x a ,y a ) = (-15,-15) and (x b ,y b ) = (30,30). Equations given in Appendix B are used in the 
PML domains for this example. Contours of p and u are shown in Figure 8 and 9 and comparisons 
with the exact solution along the line x = y are shown in Figure 10. Again, good agreement is 
found. 

To study the effects of the PML domain thickness on the reflection in computation, numerical 
solutions are computed with different PML domain thicknesses and compared with a reference 
solution. The reference solution is obtained by using a larger computational domain so that it 
is not affected by the boundary condition. In Figure 11, the maximum difference of the pressure 
between the computed and the reference solutions along the line x = 48 is shown as a function 
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of time. This measures the magnitude of the reflected wave at the out-flow boundary. PML 
domains with 6, 8, 10 and 16 grid points have been used. It is seen that reflection is reduced as 
the PML domain thickness increases. For reference purpose, results using an asymptotic solution 
based boundary condition are also shown in Figure 11. We note that the asymptotic solution based 
boundary condition would perform better if the acoustic pulse was initiated at the center of the 
interior domain and the numerical boundaries placed farther [14]. Considering that the out-going 
waves have a magnitude around 0.1, satisfactory results can be obtain with a PML domain thickness 
as small as 8 grid points. 

The third example simulates the reflection of an acoustic pulse by a solid wall with a mean 
flow M = 0.5. The computational domain is [—60,60] X [0,110] and the solid wall is located at 
y = 0, where v = 0 is applied [13]. A PML domain of 10 grid points has been used around the 
other three sides. At / = 0, an acoustic pulse of half width 5 is initiated at (x a ,y a ) = ( — 25,30). 
Since the solution can not be treated as being periodic in the y direction in this example, backward 
difference schemes of [15] are used for computing y derivatives near the upper and lower boundaries. 
Pressure contours are shown in Figure 12. It is seen that, although the pulse width is larger than 
that in the first two examples, the results are satisfactory when the same PML domain thickness is 
used. The reflection error is also assessed in Figure 13 in which the maximum pressure difference 
of the computed and the reference solutions along the line x = 48 is plotted as a function of time. 
Again, the numerical reflection is reduced as the PML domain thickness increases. The asymptotic 
solution based boundary condition, however, gives a larger reflection in this case. 

6. CONCLUSIONS 

A perfectly matched layer is proposed for absorbing the linear waves of the Euler equations. 
Wave propagation and absorption properties are given for the PML equations. It is shown that 
a PML domain so defined is capable of absorbing out-going acoustic, vorticity and entropy waves 
with no reflection (theoretically). Moreover, numerical examples that demonstrate the validity of 
the proposed PML equations are presented. 

We remark that, unlike boundary conditions based on asymptotic solutions, it is not necessary 
to apply the PML equations at far held. Although it results in solving more equations in the 
PML domains, the extra work is justified by a reduced size of the necessary computational domain. 
The PML equations also apply where an asymptotic solution is not available. In addition, the 
no-reflection and absorption properties are independent of the frequency of the waves. This implies 
that the required layer thickness may be determined independent of the wavelength of the out-going 
waves. The PML technique also suggests the possibility of using periodic conditions for computing 
spatial derivatives in certain cases, such as those in the numerical examples. 


11 



APPENDIX A : LINEAR WAVES OF THE EULER EQUATIONS 

Plane wave solutions of (1.1)-(1.4) can be expressed as follows [12] : 


( U \ 

v 

P 

\p) 


= A' 


( COS(f>\ 
sin <p 
1 

V 1 / 


g ttlJ (l+Mcos^ T 1-fMcos^ 


+ 1 + M eta t ) 


( — 5' sin ip \ 
B' cos ip 
0 

V C' ) 


g ( Aer + m“ V *) 


in which A ', B' and C represent the magnitudes of the acoustic, vorticity and entropy waves, 
respectively. 

It is, of course, possible to rewrite the above as a solution to the PML equations with a x = a y = 
0. The process is straightforward using (4.1)-(4.8) and (5.1)-(5.6). In particular, we note that the 
vorticity and entropy waves do not contribute to the n 10 and n 10 components. Moreover, by (5.1)- 
(5.2), we have u 10 /v 10 = (u 10 + u 20 )/(v 10 + v 20 ). Then, it follows that u 10 /v 10 = u acoustlc /v acoustlc . 
Thus, iiy/bo is real for solutions of the Euler equations. 


APPENDIX B : PML EQUATIONS FOR A MEAN FLOW IN GENERAL DIRECTION 

For a mean flow not parallel to the avaxis, the linearized Euler equations are 


du 

~di 

dv 

di 

dp 

di 

dp 

dt 


+ M x 
+ M x 
+ M x 
+ M x 


du 

dx 

dv 

dx 

dp 

dx 

dp 

dx 


+ M y 

+ My 
+ M y 

+ My 


du 

dy 

dv 

dy 

dp 

dy 

dp 

dy 


dp 

dx 

dp 

dy 

(du 

V<9a; 

/ du 

\dx 


dv 

dy 

dv 

dy 


(B 1 . 1 ) 

(51.2) 

(51.3) 

(51.4) 


where M x and M y are the mean velocities in the x and y directions respectively. 

It has been found necessary to split u and v into 3 sub-components. We propose the following 
PML equations : 


du 1 




d(p 

dt 

A 

O X 

U\ = 


du 2 




~M X - 

dt 

A 

O X 

u 2 = 

du 3 




~My - 

dt 

A 

Gy 

u 3 = 

dv 1 




9{p 

dt 

A 

Gy 

■Vi = 


dv 2 




-M x - 

dt 

A 

O X 

■V 2 = 


dx 

d(ui A u 2 A u%) 
dx 

d(ui A u 2 A u%) 
dy 

1 + P 2 ) 
dy 

d(v i A V’ 2 A v 3 ) 
dx 


(52.1) 

(52.2) 

(52.3) 

(52.4) 

(52.5) 
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(52.6) 


dv 3 

dt 


VyV 3 


dpi 

dt 


+ °xPl 


dp 2 

dt 

dpi 

dt 


+ °yP2 
+ ^xPl 


dp 2 

dt 


+ (J yP2 


M, atvi + : ,3+t ’ 3) 


<9(m! + u 2 + u 3 ) d(pi + p 2 ) 

dx dx 

(52.7) 

d{v 1 + v 2 + v 3 ) 9(pi+p 2 ) 

dy " y dy 

(52.8) 

d(ui + u 2 + u 3 ) d(pi + p 2 ) 

dx dx 

(52.9) 

d(vi + v 2 + v 3 ) d(pi + p 2 ) 

dy " y dy 

(52.10) 


The solutions of (B2.1)-(B2.10) can be found analogously following section 3. In particular, 
the wavenumbers k x and k y of the acoustic solution are found to be 


/ oj + iu x oj + io y 

k r = ; — — ; — - cos <p , k„ = — ■ — — ; — - sin <p 

1 + M x cos (f) + M y sin (j) 1 + M x cos f + M y sin f 

It can again be established that reflection is null at interfaces. 

We also note that when there is no mean flow, i.e. , M x = M y = 0, it is not necessary to split 
u and v. This will result in six equations for the PML domain by deleting (B2.2), (B2.3), (B2.5) 
and (B2.6). 
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Figure 1. A schematic of computational domain showing the interior domain and PML domains 
on the boundary. 
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Figure 2. Absorption coefficients on a corner of the computational domain. 
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Figure 11. Maximum pressure difference of the computed and the reference solutions along x = 48 
as a function of time. M = 0.5. Indicated are the number of grid points in PML domains used. 
Dotted line is the result when an asymptotic solution based boundary condition is used for |.t|, 
\y\ > 50. 
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Figure 13. Maximum pressure difference of the computed and the reference solutions along x = 48 
as a function of time. M = 0.5. A solid wall is located at y = 0. Indicated are the number of 
grid points in PML domains used. Dotted line is the result when an asymptotic solution based 
boundary condition is used. 
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